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qq Abstract. A discrete Fourier analysis on the fundamental domain f2<; of the d-dimensional 

lattice of type Ad is studied, where O2 is the regular hexagon and is the rhombic dodeca- 
hedron, and analogous results on ci-dimensional simplex are derived by considering invariant 

^s^j and anti-invariant elements. Our main results include Fourier analysis in trigonometric func- 

tions, interpolation and cubature formulas on these domains. In particular, a trigonometric 
Lagrange interpolation on the simplex is shown to satisfy an explicit compact formula and 
the Lebesgue constant of the interpolation is shown to be in the order of (log n) d . The basic 
trigonometric functions on the simplex can be identified with Chebyshev polynomials in sev- 
eral variables already appeared in literature. We study common zeros of these polynomials 
and show that they are nodes for a family of Gaussian cubature formulas, which provides 

1 1 only the second known example of such formulas. 
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1. Introduction 

The classical discrete Fourier analysis works with the periodic exponential functions e 2mkx 
for k in the set Z„ := {0, 1, .. . ,n — 1}. The central piece is the discrete Fourier transform 
defined via an inner product over the set of equally spaced points {k/n : k S Z„} in [0, 1), 
which agrees with the continuous Fourier transform for trigonometric polynomials of degree at 
most n. This is equivalent to Gaussian quadrature formula for trigonometric polynomials, and 
it can be used to define trigonometric interpolation based on the equally spaced points. All 
three quantities, discrete Fourier transform, quadrature, and interpolation, arc important tools 
in numerous applications and form an integrated part of the discrete Fourier analysis. 

In several variables, periodicity of functions can be defined via a lattice L, which is a discrete 
group described by AZ d , where A is a nonsingular matrix. A function / is periodic with respect 
to L = AZ d if f(x) = f{x + Ak) for any k € Z d . The usual multiple Fourier analysis of M. d 
uses periodicity defined by the lattice 1 d . One can develop a discrete Fourier analysis in several 
variables with a lattice periodicity, which has been defined in connection with signal processing 
and sampling theory (see, for example, [TUJ [14j |22]). The domain on which the analysis takes 
place is the fundamental domain of the lattice, which is a bounded set that tiles R d in the 
sense that f2 + AZ d = M. d . It is known that the family of exponentials {e 27Vla x ■. a <E L- 1 }, in 
which L 1 - is the dual lattice, forms an orthonormal basis for L 2 (f2) and these exponentials are 
periodic with respect to L. These two facts make it possible to establish a complete analogue 
of the classical Fourier analysis. Recent in [18], a discrete Fourier transform is defined and used 
to study cubature and trigonometric interpolation on the domain £1, which provides a general 
framework for the discrete Fourier analysis with lattice tiling. Detail studies are carried out 
in |18j for the hexagon lattice, for which the fundamental domain is a regular hexagon, when 
d = 2, and in [19 for the face-centered cubic lattice, for which the fundamental domain is a 
rhombic dodecahedron, when d = 3. The regular hexagon and the rhombic dodecahedron are 
invariant under the reflection group A2 and ^3, respectively. In the present paper we shall 
consider the lattice of type Ad for all d > 2, for which the fundamental domain is the union of 
the images of a regular simplex in M. d under the group action. 

Our goal is then to establish a discrete Fourier analysis on the fundamental domain of the Ad 
lattice in R d that is comparable with the classical theory in one variable as outlined in the first 
paragraph of this introduction. However, this is only the first step. The frame work developed 
in |18j applies to the fundamental domain f2 that tiles R d without overlapping, which means 
that £1 contains only part of its own boundary. For example, if the lattice is Z d , then the domain 
is the half open and half closed [— |, ^) d . The points on which the discrete Fourier transform is 
defined consist of lattice points in nfi, which is not symmetric under the reflection group Ad, as 
part of the boundary points are not included. Thus, as a second step, we will develop a discrete 
Fourier analysis that uses a symmetric set of points by working with the congruence of the 
boundary points under translation by the lattice and under the action of the group Ad, which 
requires delicate analysis. Furthermore, by restricting to functions that are invariant under the 
group, the results on the set of symmetric points in the step 2 can be transformed to results on 
a simplex that makes up the fundamental domain. This is our third step, which establishes a 
Fourier analysis on trigonometric functions on a regular simplex in ]R d . For the classical Fourier 
analysis, this step amounts to a discrete analysis on cosine and sine functions. We will define 
analogues of cosine and sine functions on a simplex in the step 3. The classical Chebyshev 
polynomials arise from cosine and sine functions. As our fourth step, we define generalized 
Chebyshev polynomials of the first and the second kind from those generalized cosine and sine 
functions, respectively, and study the discrete analysis of these algebraic polynomials. 

In each of the steps outlined above, we will develop a concrete discrete Fourier analysis 
associated with the Ad lattice in E d , providing explicit formulas and detailed analysis. The 
cases d — 2 and d = 3 studied in [TH] and provide a road map for our study. However, 
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the extension from the cases of two and three variables to the general d-variables is far from 
a trivial exercise. In order to carry out our program on the discrete analysis, it is essential 
that we have a complete understanding of the boundary of the fundamental domain and its 
congruent relations under the group. For the step 3, we also need to understand the boundary 
of the simplex and its relation with the fundamental domain. In K 2 and IR 3 , we can grasp 
the geometry of the hexagon and the rhombic dodecahedron by looking at their graphs and 
understand the relative positions between faces, edges and vertices, which are all parts of the 
boundary of the domain. This is not going to be the case in R d for d > 3. For R d we will have 
to rely on the group theory and describe the boundary elements by working with subgroups. 
The study in the cases of d = 2 and d = 3 has revealed much of what can be done and they 
will serve as examples throughout the paper. 

The generalized cosine and sine functions on the simplex, as well as the generalized Chcbyshcv 
polynomials, have been studied before in the literature. In fact, in the case of d = 2, they are 
first studied by Koornwinder in [15] , who started with the eigenfunctions of the Laplace operator 
and showed that the generalized Chebyshev polynomials are orthogonal on the region bounded 
by the Steiner's hypocycloid. Later these polynomials are generated and studied by several 
authors (see [TJ |5J |3J [7J [HI M H2J 121] and the reference therein) , often as eigenfunctions of certain 
differential operators. In [3], they appear as eigenfunctions of the radial part of the Laplacian- 
Beltrami operator on certain symmetric spaces. We highly recommend the article |3J, which also 
includes a detailed account on the history of these polynomials. Furthermore, the generalized 
Chebyshev polynomials belong to large families of special functions (see, for example, [4]), 
just as the classical Chebyshev polynomials are the special cases of Jacobi polynomials. This, 
however, appears to have little bearing with our study in this paper. It should be mentioned, 
however, that our aim is to develop a discrete Fourier analysis on the generalized Chebyshev 
polynomials, which come down to study common zeros of these polynomials, a topic that does 
not seem to have been studied systematically before. 

There are other types of discrete Fourier analysis in several variables. For the multiple 
Fourier analysis, with lattice Z d , there is a rich collection of literature as it is essentially the 
tensor product of one dimensional results. Among others, we mention the recent work on anti- 
symmetric exponential and trigonometric functions by Klimyk and Patera, see |161 117] and the 
references therein. 

The paper is organized as follows. In Section 2 we shall recall the basic facts on lattice 
and tiling, and describe the framework for the discrete Fourier analysis established in [IS]. In 
Section 3 we study the case of Ad lattice and carry out both the step 2 and the step 3 of 
our program. The best way to describe the Ad lattice and its fundamental domain appears to 
be using homogeneous coordinates, which amounts to regard M. d as the hyperplane t\ + . . . + 
td+i — of In fact, the group Ad becomes permutation group Sd+i in homogeneous 

coordinates. In Section 4 we establish a discrete Fourier analysis on the simplex, which studies 
trigonometric polynomials on the standard simplex in R d . One of the concrete result is a 
Lagrange interpolation by trigonometric polynomial on equal spaced points on the simplex, 
which can be computed by a compact formula and has an operator norm, called the Lebesgue 
constant, in the order of (logn) d . Finally, in Section 5, we derive the basic properties of the 
generalized Chebyshev polynomials of the first and the second kind, and study their common 
zeros, which arc related intrinsically to the cubature formula for algebraic polynomials. In 
particular, it is known that a Gaussian cubature exists if and only if the set of its nodes is 
exactly the variety of the ideal generated by the orthogonal polynomials of degree n, which 
however happens rarely and only one single family of examples is known in the literature. We 
will show that this is the case for the generalized Chebyshev polynomials of the second kind, 
so that a Gaussian cubature formula exists for the weight function with respect to which the 
Chebyshev polynomials of the second kind are orthogonal, which provides the second family of 
examples for such a cubature. 
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2. Discrete Fourier analysis with lattice 

In this section we recall results on discrete Fourier analysis associated with lattice. Back- 
ground and the content of the first subsection can be found in [6j [10l HH [22] . Main results, 
stated in the second subsection, are developed in [18J but will be reformulated somewhat, for 
which explanation will be given. We will be brief and refer the proof and further discussions to 



2.1. Lattice and Fourier series. We will consider lattices on M. d but will use homogeneous 
coordinates, which amounts to identify ~R d with the hyperplane X\ + . . . + x d+ \ = in R d+1 . 
For convenience, we will state the results on the lattices on a d-dimensional subspace V d of the 
Euclidean space K m with m> d. Since V d can be identified with K d , all results proved for M. d 
can be restated for V d . 

Let Y d be a ^-dimensional subspace of M. m with to > d. A lattice L of Y d is a discrete 
subgroup that contains d linearly independent vectors, 

L := {kiai + k 2 a 2 + • • • + k d a d : ki G Z, i = 1, 2, • • • , d} , 

where a\, ■ ■ ■ ,a d are linearly independent column vectors in Y d . Let A be the matrix of size 
to x d, which has column vectors a±, ■ ■ ■ ,a d . Then A is called a generator matrix of the lattice 
L. We can write L as AZ d ; that is 

L = AZ d := {Ak : ke Z d } . 

The dual lattice L 1 - of L is given by 

L 1 - := {x e Y d : x ■ y e Z for all y £ L} , 

where x ■ y denotes the usual Euclidean inner product of x and y. The generator matrix of 
L- 1 , denoted by A , is given by A 1 - = A(A tr A)^ 1 , which is the transpose of the Moor-Penrose 
inverse of A. In particular, if m = d, then the generator matrix of L is simply A~ tr and 
L 1 - = A- tr Z d = {A- tr k : k£Z d }. 

A bounded set c Y d is said to tile V d with the lattice L if 

Xn(x + a) = 1 for almost all x £ V d , 

where xn denotes the characteristic function of f2, which we write as + L = Y d . Tiling and 
Fourier analysis are closely related as demonstrated by the Fuglede theorem. Let J n f(x)dx 
denote the integration of the function / over il. Let (•, -)q denote the inner product in L 2 (il), 

(2-1) (f,g) a ~±Jf(x)tfx)dx, 

where |0| denotes the measure of f2 and the bar denotes the complex conjugation. The following 
fundamental result was proved by Fuglede in [13J. 

Theorem 2.1. Let c V d be a bounded domain and L be a lattice ofY d . Then fi + L ~Y d if 



and only if |e 27r4K ;E : k £ L^} is an orthonormal basis with respect to the inner product (2.1 1 
Written explicitly, the orthonormal property states that 

(2.2) _L^e 2 "^dx = <5 K , , ^i 1 



Moreover, if L = AZ d , then the measure |0| = v /det(^ tr A). 

The set is called a spectral set for the lattice L. If L = AZ d we also write Ha instead of £1. 
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Because of Theorem |2.l| a function / 6 can be expanded into a Fourier series 

The Fourier transform / of a function defined on L 1 (N d ) and its inversion are defined by 

/(£) := / f(x)e ™*<°dx, /(*) := / /(Oe 2 ™^^. 
Our first result is the following sampling theorem (see, for example, [HI |2"2"| ). 

Proposition 2.2. Let O be the spectral set of the lattice L. Assume that f is supported on 
and / G L 2 (0). Then 

/(*)= ^ /(«)*„(*-«) 

in I/ 2 (0), where 

This theorem is a consequence of the Poisson summation formula. We notice that 

$fi(«0 = <5o,k for all keL 1 

by Theorem |2.1[ so that 4>n can be considered as a cardinal interpolation function. 

2.2. Discrete Fourier analysis and interpolation. A function / defined on V d is called 
periodic with respect to the lattice L if 

f(x + a) — f(x) for all a E L. 

The spectral set O of the lattice L is not unique. In order to carry out the discrete Fourier 
analysis with respect to the lattice, we shall fix an O such that O contains in its interior and 
we further require that O tiles V d with L without overlapping and without gap. In other words, 
we require that 

(2.3) ^xn(£ + a)=l for all x G Y d . 

For example, for the standard cubic lattice Z d we can take O = [— |, |) . Correspondingly a 
spectral set satisfying (2.3 1 is called a fundamental domain of Lattice L. 

Definition 2.3. Let La = AZ d and Lb = LVL d be two lattices ofV d , £Ia and O^ satisfy (2.3 1. 
Assume N := B 1 ' A is a nonsingular matrix that has all entries being integers. Define 

A N :=QA<^Lg and Ajy- := Ob n 

Two points i,i/6 V d are said to be congruent with respect to the lattice AZ d , if x — y G AZ d , 
and we write x = y (mod A) or x = y (mod AZ d ). The following two theorems are the central 
results for the discrete Fourier transform. 



Theorem 2.4. Let A, B and N be as in Definition 2.3 Then for any k G L\, 

1 V^ r 2« K .a_/ 1 > if k = (mod B), 



^'^ | det (AT) | ^ 6 If), otherwise, 



and 



(2.5) 



1, if a = (mod A), 



| det(A r ) | | 0, otherwise. 
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Theorem 2.5. Let A, B and N be as in Definition \2.3\ Define the discrete inner product 

{L9)N = \d^m £ f{a)W) 

1 y n aeAjv 

for f, g € C(VLa), the space of continuous functions on VLa- Then 
(2-6) </, g) N = (/, 5 ) 

/or all f, g in the finite dimensional subspace 

U N := span {<f> K : <f> K {x) = e 2 ™ K ' x , KeA^j. 



Let denote the cardinality of the set E. Setting k = or a = in ( 2.4 1 or (2.51, 
respectively, we see that 

(2.7) lAjvlHAjvtMdetCAT)!. 

In particular, the dimension of Hn is |Ajy| = | det(iV)|. 

Let Xjv/ denote the Fourier expansion of / £ (7(0^) in 7~Ln with respect to the inner product 
(•, •) m- Then, analogous to the sampling theorem in Proposition 2.2 2jv/ satisfies the following 
formula 

Infix) = J2 /(«)*n B («-«). /e^a), 

aGA N 

where 

«*> = RsW E 

The following theorem shows that Iat/ is an interpolation function. 



Theorem 2.6. Lef A, -B and N be as in Definition \2.S\ ThenXj^f is the unique interpolation 
operator on N in TLn, that is 

Z)v/(a) = /(a), Va G Ajy- 
In particular, |Ajv| = |A]y|. Furthermore, the fundamental interpolation function 3>$f B satisfies 

(2.8) *n s ( a; )= E $ n fl (* + *)• 

kEL/1 



Proof. Equation ( 2.8 1 was proved in [TH] . Using ( 2.5 1 gives immediately \&q b = <5 a .o for a S , 
so that the interpolation holds. This also shows that {^n — a) : a £ Ajy} is linearly 

independent. By ( |2~7| , |Ajv| = |det(iV)| = \f^ N \. Furthermore, ( |2~4| and ( |2~5| show that the 
interpolation matrix M = (4> K (ct)) a£ j^ N is invertible. Consequently the interpolation on 
points in Ajv is unique. □ 

This theorem is stated in |18j under the additional requirement that A jy = A^r" , which holds 
in particular in the case that A is a constant multiple of B. We prove the more general version 
here for future references. 

The results state in this subsection provide a framework for the discrete Fourier analysis 
on the spectral set of a lattice. We will apply it to the d-dimensional lattice associated with 
the group Ad in the following section. As mentioned in the introduction, the case d — 2 and 
d = 3 have been considered in detail in [TH] and [IS], respectively. These lower dimensional 
cases provide a roadmap for the d-dimcnsional results. 
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3. Discrete Fourier analysis on the fundamental domain 
3.1. Lattice A d and Fourier analysis. For d > 1, we identify R d with the hyperplane 

< +1 := {{tx,t 2 ,--- ,t d+1 )eR d+1 :t 1 +t 2 + ---+t d+1 =0} 
of R d+1 . The lattice A d that we will consider in this paper is simply constructed by 

Z^ 1 := Z d+1 n M^ 1 = {(ki,k 2 , ■ ■ ■ ,k d+1 )eZ d+1 :k 1 + k 2 + --- + k d+1 = Q}. 

In other words, we will use homogeneous coordinates t G R 1 ^ 1 to describe our results in d- 
variables and A d — Zf^ 1 in our homogeneous coordinates. 

Throughout this paper, we adopt the convention of using bold letters, such as t and k, to 
denote homogeneous coordinates. The advantage of homogeneous coordinates lies in preserva- 
tion of symmetry. In fact, many of our formulas are symmetric and more transparent under 
homogeneous coordinates [18j HU |27j [28] . 

The lattice A d is the root lattice of the reflection group A d [HI Chapter 4] . Under homoge- 
neous coordinates, the group A d is generated by the reflections {<7jj :l<i<i<rf+l}, where 
Oij is defined by 

toy := t - 2- — ' . e itj = t - (U - t^e^j, where e hJ := e, - ej. 

The last equation shows that cx^ is a transposition. Thus, the group A d is exactly the permu- 
tation group Sd+i of d+ 1 elements. Let := S<j+i- For t S M^ 1 and cr G the action of rr 
on t is denoted by tcr, which means permutation of the element in t by a. 

The fundamental domain that tiles Wj^ 1 with lattice ^ is the spectral set of A d , 

(3.1) n H := {t G : -1 < U -tj < 1, 1 < i <j < d + l}. 

The strict inequality in the definition of £Ih reflects our requirement in ( |2.3[ ) . In the case of d = 2 
and d = 3, the fundamental domain are hexagon and rhombic dodecahedron, respectively, and 
they are depicted in Figures |3.1||3.2| in which vertices are labeled in homogeneous coordinates. 

(-1. 1,1)/ aHH' D 



V 3' 3' 3 



FIGURE 3.1. Fundamental domain — the regular hexagon for the lattice A 2 . 

As described in [HJ Chapter 21], the spectral set is the union of its fundamental simplex 
Ah, defined by 

(3.2) A H := {teK^ 1 : < t { - tj < I, l<i<j<d+l}. 



8 



HUIYUAN LI AND YUAN XU 




Figure 3.2. Fundamental domain — the rhombic dodecahedron for the lattice A 3 . 



Let t, t, . . . , t be denoted by {t} . The d non-zero vertices of A H are given by 

k 

^-y, l<k<d, where v fc := ({d + I - k} k , {-k} d+1 - k ) . 

Let tQ denote the orbit of t under Q, that is, tQ := {tcr : <r E Q}. Then is the union of the 
images of A H under Q, that is 

(3.3) Q H = A H Q:= \J tg= \J {ta:te A H } . 

teA H aeg 

Furthermore, the partition is non-overlapping, i.e., for any t,s € Ah and t =/= s, tQ n stj = 0. 

For any t £ M^ 1-1 , the stabilizer of t is denoted by Qt '■= {<? £ Q : to = t}. Since Q is finite, 
\Qt\ x |t5| = \G\. For the vertices of A H , \v k G\ = k $£}! k) , and |0 V *| = fc!(d+ 1 - fe)!. It is 
worthwhile to note that fljj is the convex hull of the Y^t=i W k G\ = — 2 vertices in 

{^■■°^,l<k<d}. 

The generator matrix of the lattice Ad is given by the (d + 1) x d matrix 



A := 



/ 1 


. . 


. 


°\ 





1 . . 


. 








. . 


1 








. . 


. 


1 




-1 . . 


. -1 


-v 
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and the generator matrix of the dual lattice Lj[ is the (d + 1) x d matrix 



A 1 - = AiA^A)' 1 



/ d 


-1 . . 


. -1 


-1\ 


-1 


d . . 


. -1 


-1 


-1 


-1 . . 


. d 


-1 


-1 


-1 . . 


. -1 


d 


V-i 


-1 . . 


. -1 


-v 



To describe explicitly, we introduce the notation 

H := {k G Z^ +1 : ki = k 2 = ■ ■ ■ = k d+1 (mod d + 1)} . 



For j G Z set k = (d+l)A j. It follows readily that k G Z^ and a quick computation shows 
kel. On the other kind, the definition of A- 1 leads to j = A tr k/ (d+ 1), which shows j G Z 3 
if k € EL Consequently, the dual lattice is given by 

k 



1 



k g 



For e 2 ™ 1 ^* with k G L±, we introduce the new notation 



(3.4) 



e d + 



U-* for j e 



By the Fuglede theorem, {0j : j G H} forms an orthonormal basis in L 2 (f2#); that is, 



(3.5) 



&k,</>j) = 5kj, k,j e 



where the inner product is defined by 

1 



(3.6) 



h I Jn 



/(t)fl(t)dt 



1 



/(t). 9 (t)rft. 



Vd+1 

Definition 3.1. A function f is H-periodic if it is periodic with respect to the lattice A^; that 
is, f(t) = /(t + k) for all k G Z^ +1 . 



1 



/(t)0_ k (t)dt. 



Evidently, the functions 0j(t) in (3.4 1 are ff-periodic. Furthermore, (3.51 shows that an 
iJ-periodic function / can be expanded into a Fourier series 

(3.7) J~^jy>k(t), where / k := — L 

keH 

Recall that s = t (mod H) means s — t G Z 1 ^ 1 . The following lemma will be useful later. 
Lemma 3.2. J/t,s G SIh and s = t (mod H), then t = s. 

Proof. If t, s G and t = s (mod H), then s-te Z^ 1 " 1 and, set k := s — t, —1 < ki — kj < 1 
for all 1 < i,j < d + 1. The last condition means that either ki G {0, 1} for all 1 < i < d+ 1 
or fcj G {0, —1} for all 1 < i < d + 1. The homogenous condition = then shows that 

k = or s = t. □ 

3.2. Structure of the fundamental domain. In this subsection, we concentrate on the struc- 
ture of Hh- For this purpose, we set N^+i := {1, 2, • • • , d+ 1} and start with the observation 
that CIh can be partitioned into d + 1 congruent parts. 

Lemma 3.3. For 1 < j < d + 1 define 

:= {t G : < U - t } < 1 and < - tj < 1 /or 1 < i < j < Z < d + l} . 

Then 



and 



n { ^ } n = /or 1 < i ^ j < d + 1. 



l<j<d+l 
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Proof. For t G f2#, let J = {k £ N^+i : tk < U for i G N^+i} and j = mitifc 6 j k. The definition 
shows that ti — tj > for all i G N^+i; furthermore, if i < j < k then, as t G £Ih , < ti — tj < 1 



and —1 < £j — < 0, so that t G which implies that C UjeN d+ i ^ff' ■ Since 

fijj* C f2# for each j G N<j+i, it follows that Cl H = (J . 



If t G f2^ } n with 1 < % < j < d + 1, then t G so that < ^ ■ - ^ and t G so 



that < ti — tj. As a consequence, tj < ti < tj, which implies D Q,^p 



□ 



3.1 



the sets Q { ^ } ,Q { ^ } ,Q { ^ } 
are the three parallelograms inside the hexagon, starting from the one in the upper left rotating 



For d = 2, the partition of the hexagon is evident from Figure 



clockwise. In the case d = 3, the rhombic dodecahedron is more complicated (see Figure 3.2 1, 



its decomposition is depicted in the Figure 3.3 



(fi-4.~l.-i) 



(0,0,0,0) 





(-1 -1 -1 5) 

V 4 ' 4 ' 4 ' 4 > 




Figure 3.3. Rhombohedral partition of the rhombic dodecahedron. Upper- 
Left: f4 1} ; Upper-Right: f^ 2} ; Lower-Left: Q^; Lower-Right: Q.^ . 

The next lemma clarifies the intersection of the closure of Qn\ for which we need to define 
il H for a set J C N^ + i. 
Lemma 3.4. For C J C Nd+i, define 

n J H := {t G K# +1 : U = tj, G J; and 0<U-tj< L Vj G J, Vi G N d+ i \ j} . 

T/ien 

iln = and fl H = &]P ■ 

Proof. The first equation follows easily from the previous lemma by taking closure. 



If t G n then < ti — tj < 0; that is, ^ = tj. Hence, if t G Plje/^fl"' then 

ti = tj,yi,j G J, which implies Hie/^^ £ ^ff- Since f2^- C fijij^ whenever j G J by 



do} 



{/} 



definition, we conclude that 



0} 



□ 
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3.3. Boundary of the fundamental domain. In order to carry out the discrete Fourier 
analysis on the fundamental domain fin, we need to have a detailed knowledge of the boundary 
of the region. Much of this subsection is parallel to the study in [T5]> where the case d = 3 is 
discussed in detail. Some of the proofs, in fact, are essentially the same as in the case of d = 3 
and often only minor adjustment is needed. In such cases, we shall point out the necessary 
adjustment and refer the proof to [H?] . 

We use the standard set theoretic notations dtt, tl° and to denote the boundary, the 
interior and the closure of f2, respectively. Clearly ft = Q° U dfl. For i, j £ N^+i and i =^ j, 
define 

Ftj = {t£n H :U- tj = 1}. 

There are a total d(d+ 1) distinct Fij, each stands for one facet of d— 1 dimension, together 
with its boundary, of the fundamental domain. 

The boundary of our d-dimcnsional polytope consists of lower dimensional sets, which can 
be obtained from the intersections of d — 1 dimensional facets. For example, for d — 3, the 
boundary consists of faces, edges, and vertices; edges are intersections of faces and vertices are 
intersection of edges. In order to describe the intersections of facets, we define, for nonempty 
subsets /, J of Nd+i, 

Qjj ■= P| ={te U H : tj =t l -\ for all i £ I,j £ J} . 

The main properties of these sets are summarized in the following lemma, the proof of which 
is given in [15] . 

Lemma 3.5. Let I , J, Jj, Jj be nonempty subsets o/N^+i- Then 

(i) Q ItJ = if and only if I D J ^ 0. 

(ii) QxuJx n Qi 2t j 2 = if hUl 2 = I and J 1 (JJ 2 = J. 

To carry out a discrete Fourier analysis, we need to distinguish points on the boundary 
elements of different dimensions. It is also necessary to distinguish a closed boundary element 
and an open one. For example, for d = 3, we will distinguish a face with its boundary edges 
and a face without its edges. To make these more precise, we introduce the notation 

lC = {(I,J):I,JcN d+1 ; in J = <D}, 

/C = {(I, J) £ K. : i < j, for all (i,j) £ I x J} . 

Definition 3.6. For (I, J) £ K,, the boundary element Bi j of the fundamental domain, 

B ItJ := {t G n T>J : t £ Q htJl for all (h, J x ) £ K with \I\ + |J| <\h\ + \h\} , 

is called a face of dimension k, or simply a k-face if k = d+ 1 — |/| — | J|. In particular, it is 
called a facet if k = d — 1 (or \I\ + \ J\ = 2), a ridge if k = d — 2 (|/| + \ J\ — 3), a face if k = 2 
(\I\ + |J| =d-l), an edge if k = 1 (|/| + |J| = d) and a vertex if k = (|/| + | J\ = 

Figure [5~4] shows a 3-face, or a facet, of the fundamental domain when d = 4. 

In the following, when we say a fc-face we mean the open set, that is, without any of its 
boundaries. For fc-faces with k > 1, the boundary elements Bj j represent the interiors. In 
fact, it is easy to see that B^ ^ = F?- for distinct integers i, j £ Nd+i and Bj^j = 7 for 
(/, J) £ K. By the definition of Q,h and we can write Bi j more explicitly as 

(3.8) B hJ = {t£ : U>ti> tj =t t -l for all i £ I, j £ J and I £ N d+1 \ (I U J)} . 
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Figure 3.4. The 3-face (facet) t\ — t 5 = 1 of the fundamental domain when d = 4. 



The boundary elements of the same type will often be considered as a group. For this 
purpose, we further define, for < i, j < i + j < d+ 1, 



(3.9) 



K**:={(I,J)€1C: \I\ = i, \J\=j], l'> ' : 
{(I,J)eJC : \I\=i, \J\=j], Bl' j 



1(1,3 . 



U B ItJ 

= U B Jt 

(/,,/)e/cj, 3 



Proposition 3.7. Let (I, J ) E K. and (I\,Ji) E JC. 
(i) B h j n B IuJl = 0, ifl^h and J ^ J x . 



(ii) n ff \^ = u (J n eK v = U 



- l 0<ij<i+i<d+l 



(in) n H \ n° H - U(/,j) 6 Ko Sj v _ Uo<i,j<i+i<d+i 

The proof of this proposition is essentially the same as the one given in [13] for d — 3. 

The boundary elements can also be described by using symmetry. In fact, boundary elements 
of the same type can be transformed by the group Q = Sd+i- To make it more precise, we define, 
for (I, J) 6 K and a E Q, 

Bi,j<t := {ttr : te%}. 

Then it follows readily that 

B \i\,\J\ = (J flj J(T := {ta:te Blj, ueG}. 
teg 

Recall that er.y = t — (t , — tj)e it j denotes the transposition in C/ that interchanges i and j. 
We clearly have er^ = cr^j and ajj is the identity element of the group. For a nonempty set 
I C Nd+i, define Qj := {aij : i,j E I}. It is easy to verify that Qi forms a subgroup of Q = Sd+i 
of order |/|. 

Lemma 3.8. Letfj^B 1,3 and#B l J denote the number of distinct (d+1 — i — j)- faces contained 
in B l ' J and B^ 3 , respectively. Then 

(d+iy. (d+iy. 



i\j\(d+l-i-j)\ 



and ^Bq 3 



(i + j)\(d+l-i-j)V 
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Proof. By definition, _B lJ and Bq 3 are unions of (d+ 1 — i — j)-faces of the fundamental domain. 
The first formula follows from the fact that Bjja = Bi j if a G Qj U Qj U 5n d+1 \(/uJ)- The 
second one follows from the fact that , for any K C Nd+i and i, j > 1 such that \K\ = i + j, 
there is a unique (I, J) G /Cg J such that I U J = K. □ 

Using these formulas, we can derive by setting k = d+X — i — j that the fundamental domain 
£Ih has a total 



V _ r2 d+i-fc o^ rf+ M 0<fc<d-l 

Ajj!Jfc!(d+l-i-*)!~^ } { k )' °< fc < d 1. 

distinct k— faces, while f2# has a total (d— fc)^^ 1 ) distinct k— faces. 

For the periodic functions, we will need to consider points on the boundary elements that 
are congruent modulus Z^ 1 . For (/, J) C K, we further define 

[Bi,j] := {flj.j + k : k G Z^ +1 } n H H = {t + k G H H : t € Bj.j.ke . 



By lemma 3.2 the set [Bij] consists of exactly those boundary elements that can be obtained 
from Bi j by congruent modulus Z^ 1 . As an example, in the case of d = 3, we have 

(3 1Q) B {i},{2,3} = {(t,t - l,t - 1,2 - 3t) : \ < t < f } , 

B { i,2},{3} = {(l-t,l-t.-*,3t-2): \<t<l), 
and from this explicit description we can deduce, for example, 

[£{1},{2,3}] = £{1},{2,3} U (£{1},{2,3} + (-1, 1, 0, 0)) U (S {1} , {2i3} + (-1,0, 1, 0)) 
= - B {1},{2,3} U £{2},{1,3} U £{3},{1,2}- 

The last equation indicates that [Bj is a union of Bji ji, which is stated and proved in [TO] 
for d = 3 and the proof can be adopted with obvious replacement of d — 3 by general d. 

Lemma 3.9. Let (I, J) G K. Then 

(3-11) [Bi,j}= (J Bj,ja. 

Since can be obtained from /C ' 3 from the action of G, it follows that 
(3.12) &*= |J [Bj,j]= |J [B], 0<i,j <i+j <d+l. 

We also note that [B h j \ n [S/^jJ = if (I, J) ^ (Ji, Ji) for (I, J) G /C and (Ji, Ji) G £ , 
which shows that (3.121 is a non-overlapping partition. 

To illustrate the above partitions, let us consider the case d = 2, for which 

B 1 ' 1 = [-B{i},{2}] U [-B {1 } !{3 }] U [-B{ 2 },{3}] 5 

B 1 ' 2 = [£{!}, {2,3}L B 2 ' 1 = [-B{1,2},{3}], 

where 

S {i},{2} = {(1 -*,2i- : | < t < §} , 

"^{1}.{3} = B{1} : { 2 }0'2,3, #{2},{3} = £{1},{2}°'1,2) 
S {1},{2,3} = {(fj-li-s)}) B {1,2},{3} = {(g, !>-§)} ■ 

In the case <i = 3, the boundary elements are given explicitly in [19 . 
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3.4. Point sets for discrete Fourier analysis for the lattice Ad- We apply the general 
result on discrete Fourier analysis in Section 2 to the lattice Ad by choosing A — A and B = nA 
with n being a positive integer. Let 1 = Then the matrix 



N 



B tr A = nl + nt tr l 



(In 


n 


n 


n\ 


n 


2n . 


n 


n 


n 


n 


. 2n 


n 


\n 


n 


n 


2n) 



has integer entries. Since N is now a symmetric matrix, An 
that AL = Ajy. We denote A^r by H„, which is given by 



An"; moreover, it is easy to see 



k 



-n < 



1 <n 1 l<i<j<d+l 



d+1 

The finite dimensional space TLm of exponentials in Theorem |2.5| becomes 

H n := span {fa: ke H n } with dimW n = det(iV) = (d+ l)n d . 

The points in H n are not symmetric under Q , since points on half of the boundary are not 
included. We further define the symmetric extension of H n by 

kj ' k o 



:= ke 



(d+l)n 



e n H = k e 



-n < 



1 



<n, 1 <i < j < d+1 



For d — 2 and 3, the scope of H* is depicted in Figure 3.5 and Figure 3.6 in which the vertices 
are labeled in homogeneous coordinates. 




Figure 3.5. H* for d = 2. 

For discrete Fourier analysis, it is essential to understand the sets H n and H* , for which the 
main task lies in studying the structure of points on the boundary; in other words, we need to 
understand exactly when ^ d +i)n belongs to a specific boundary elements of fin- Much of the 
work in the previous subsection is a preparation for this task. In the rest of this subsection, we 
develop results in this regard. 

We start with the partition of H n and H* into d+1 congruent parts, each within a paral- 
lelepiped, corresponding to the partition of and flu , as shown in Figures 3.3 for d = 3. 

Lemma 3.10. For 1 < j < define 

H#> = (k S H : 1< k \~ k J < n andO< k lZ-k < n - 1 for l< i < j < l< d + l\ . 
{ ~ d+1 ~ ~ d+1 ~ J 
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(2n,2n,-2n,-2n) 




(-2n,-2n,2n,2n) 



Figure 3.6. Ml for d = 3. 



Then 

H„ = |J H^' } and n H^' } = /or 1 < i ^ j < d + 1. 



Proof. By definition, Hn"^ = |k G H : € ^h^}- Hence, this lemma is an immediate 

consequence of Lemma |3.3| □ 



Similarly, as a consequence of Lemma |3.4| we have the following lemma. 
Lemma 3.11. For C J C N d+ i, define 

W,[ :={keH:i, = kj, Wi,j € J; and < - % < (d+l)n, Vj G J, Vi € N d+i \ J} . 

h:= U h » } anrf H « = fl H » } - 

Next we consider further partitions of H* . The set of interior points is defined by 



The set of points on the boundary of fin is then H* \ H° . Recall (3.9 1, we define, for < i, j < 
i + j<d+l, 

(3.13) E# := {k G H : ^ G B 1 ' 3 j , <f := {k G H : ^ G ///,<}. 

The index set consists of those points j in H ra such that ^ d ^ n are in the boundary element 
of dn H . Using Proposition pL7l it is easy to see that n EI^ 1 Jl = if i ^ i x or j ^ ji, 

(3.14) |J H«=H:\H° and |J H« = H n \H°. 

0<ij<i+j<d+l 0<ij<i+j<d+l 

Lemma 3.12. For n. > 1, < i, j < i + j < d + 1, 

|H°| = n d+1 - (n-l) d+ \ |H«| = J^fn-ir-'^. 

1 - j)\ 
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Proof. The first equation follows from |H° 

(d+iy. 



3.8 



B l,J consists of 



|H*_]J and (3.191 below. As shown in Lemma 
distinct boundary elements Bij with (/, J) 6 K., \I\ = i and 



\J\ = j, which are (d + 1 - 
points in H* as seen by (3.- 



j)-faces. Each of these faces contains (n — l) d+1 1 J distinct 
from which the second equation follows. □ 

We will not need to define the point set that corresponds to B I J for each pair (/, J) C /C 4,J . 
However, we do need to understand the congruent relation modulus Z^ 1 among points in H* . 
For each j € 



I^ J , we define 



5, := {k € 



k - J 

(d+ 



(3.15) ^ n. . (d+1) „ 

Since, by definition, [-Bj,j] = {s : s = t (mod Z^ 1 )^ € -£>/.,/}, it follows from (3.121 that 

H#= (J ^ 

By pTT] ), [B/,./] is the union of ( |J| ,t |J| ) components of B UJ a. It follows that = ( k f) for 



3.5. Ad Fourier partial sum. For the Fourier expansion (3.7 1 of an ii-periodic function, we 
define its partial sum with respect to the fundamental domain of Ad as 

(3.16) S n /(t) := V (/, 0k>#k(t) = I /(s)Af (t - B)ds, 

keH » Vd+lJn H 

where D„ is the Dirichlet kernel for the Fourier partial sum 

(3.17) ^(t):= £ c^". 

keH* 

The following theorem gives a compact formula for the Dirichlet kernel. 
Theorem 3.13. For n>0, 

d+l . 

(3.18) D«(t)=9 n+1 (t)-e n (t), wftere 6„(t) := TT ^ZLz.. 

- LJ - sin7rt,- 
j=i J 

This theorem is stated and proved in |19j for d = 3 and the proof carries over to general d 
with the obvious replacement of d + 1 = 4 by general d+l. Recently a different derivation of 
this formula is given in [30 . 

As one consequence of the compact expression in (3.181, we obtain 

(3.19) \W* n \=DX(0) = (n+l) d+1 -n d+1 . 

Another result that follows from the compact formula of the Dirichlet kernel is an upper bound 
for the Lebesgue constant, which is the norm of the partial sum S n f in (3.16). Let |j/||oo denote 
the uniform norm of / <E C{Q,h) and let ||SVi||oo denote the operator norm of S n : C(Q,h) i— > 
C(Q„). 

Theorem 3.14. There are positive constants c and C independent of f and n such that 

(3.20) c(logn) d < ||5 n |U < C(logn) d . 

The upper bound ||i5 n ||oo < C(logn) d was proved in for d = 3, the proof extends to the 
general d. It turns out, however, the upper bound can be derived from a general result in |24j . 
which establishes the inequality 



(3.21) 



k: k/n£E 



dx < C(logn) d 



DISCRETE FOURIER ANALYSIS IN d- VARIABLES 



17 



for E being a polyhedron in M. d . Indeed, choosing a constant a, if necessary, the set aO can 
be imbedded inside [— 77,77] d and a can be absorbed into the index set, so that ( |3.20[ ) follows 
from the general result. The lower bound that matches (3.211 was established in |31] for E 
being a convex polyhedron in W l . Upong choosing an appropriate dilation constant a so that 
ail contains [— 7r,7r] d , the lower bound in (3.201 follows. 

3.6. Discrete Fourier analysis on the fundamental domain. For the lattice Ad, the 
general result on the discrete inner product, Theorem 2.5 gives the following proposition in 
homogeneous coordinates: 

Proposition 3.15. For n > 0, define 

1 



(/, 9)n ■= 



(d+l)n d 



The 



if, g) = (/j g) n , f, g e u n . 

The inner product (■,•)„ is defined as a sum over the index set H„, which is not symmetric 
over £Ih as only part of the points in the boundary of £Ih are accounted for. More interesting 
to us is to consider an inner product based on the symmetric index set H* , which turns out to 
be equivalent to (•,•)« for iJ-periodic functions. 

Definition 3.16. For n > define the symmetric discrete inner product 



(/: g)n 



(d + 1) 



j6H* 



where c\ =1 if j S H° and c[ 



(n) 



i/j G 



"j - v j ■ n> "j 

For instance, if d = 2 then is a regular hexagon and we have 

1, j € H° , (n 3 — (ri — l) 3 points in the interior), 

j G H*' 1 , (6 x (n — 1) points on the edges), 

jeii' 2 Uiy, (6 vertices); 
if d = 3, then ri// is the rhombic dodecahedron and we have 



(n) _ 



(n) , 
cj = < 



1, 

1 

2 ' 
1 

3 ' 
1 

4 ' 
1 

k 6' 



j e 
j g 
j g 

j G 
j G 



I 1 - 1 

I 1 ' 2 

I 1 ' 3 
I 2 . 2 

"71 ' 



UH 2 ' 1 , 
UH 3 .\ 



(n 4 — (n — l) 4 points in the interior), 
(12 (n — l) 2 points on the faces), 
(2 x 12(n — 1) points on the edges), 
(2x4 points on the vertices), 
(6 points on the vertices). 



jgH*\H° 



E 



It is easy to verify that 

0<i,j<i+j<d+ 

/d+1 
V + 3 

'd 



E 

0<i,j<i+j<d+l 

(«*+i)E(^)(»-i) 
fe=i 



(n-1) 



\d-fc 



d+1 
ft=2 

d+1 

d+l-fc 



d+1 



(n-1) 



d+l-fc 



fc=2 



d+1 



(n-l)« 



(d + 1) (n d - (n - l) d ) - (n d+1 - (n - l) d+1 - (d + l)(n - If) 
(d + l)n d -n d+1 + (n-l) d+1 , 
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so that (lA)l= xs ^(\W n \ + (d+l)n d -n d+1 +(n-l) d + 1 )=lby < p9 l. 
Theorem 3.17. For n > 0, 

(f, g) = (f, g)n = </, g)* n , f,ge H n . 

Proof. Let / be an if -periodic function. Then 

».-' « ^ ^ 1 



X! C j ^( (d+l)n ) ^ X] 
jSH* \H° 



0<i,fc<i+fc<d+l 



E 



i 

i+k\ 



E /( 



(i+k\ 

0<i,k<i+k<d+l V i / jgjji.fc keSj 



E E f((d+l)n)' 



Since |<Sj| = ) for j € , using the invariance of /, we then conclude that 



E C j /((d+l)n) _ E 



jGH*\H° 0<i,fc<i+fc<t/+l 
» _ 



E 



i + fc 



Since c. = 1 if j G H° , the proof is completed. 



□ 



Since Theorem 3.17 shows that the integral of fg G 7i n agrees with the discrete sum over 
H*, it is not surprising that we have a cubature formula, which turns out to have a high 
degree of precision. To be more precise, we define by T n the space of generalized trigonometric 
polynomials, 

T n := span {0 k : k G H* } . 



Theorem 3.18. For n > 0, the cubature formula 



(3.22) 



1 



Vd + 1 
is exact for all f G T% n —\. 



/(t)dt 



1 



(d + l)n d . 



E C j f((d+l)n) 



j6H« 



Proof. It suffices to prove that (3.22) is exact for /(t) = </>k(t) with any k G E^ n _i. Since 
f2# tiles Mj 1 " 1 with the lattice A d , there exit s G and 1 G Zp' 1 such that pjrc = s + 1. 

Noting that k G and 1 G 1 are vectors with integer entries, we further deduce that 

for certain m G H n , which states that k = m + [d + l)nl with certain m G H„ and 



s = 



(d+i) 



1 G Z'j+K Next we show that m = only if k = 0. In fact k = (d + l)nl G Hjj^-i clearly states 
that G fi n- Then the restriction - - 2 < h - L < 2 - - for 1 < i < j < d + 1 and the 

2n— 1 n n 1 J — n — J — 

homogeneity of 1 imply that 1 = 0, which gives k = in return. Now applying the periodicity 
of </>j, Theorem 3.17| and (3.5 1 yields that 



(d 



+ lW E 4 Vk( (rf + 1)n ) - ( d + i) n d E C j Vm( (d + 1)n ) 



i) <Ao) — ^m,0 — <^k,0 — 



t (t)dt. 



This completes the proof. 



□ 



The cubature (3.22) is an analogue of the classical quadrature formula for trigonometric 
polynomials of one variable (see, for example, [MJ Vol. II, p. 8]), which holds for trigonometric 
polynomials of one variable and has vast applications. 
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3.7. Interpolation on the fundamental domain. For the fundamental domain of Ad, the 
general result on the interpolation, Theorem 2.6 becomes the following: 

Proposition 3.19. For n > 0, define 

2n/(t) := - Tsfei)' Where ®n(t) = , ^ d ]T 

jeH„ ^ ' keH„ 

/or / £ C(il# ). TTien Z™/ € 7i n and 

The function l n f interpolates on points in H„, which is not symmetric as noted before. We 
are more interested in interpolation on all points in H* . The operator X* that we are able to 
define, however, does not interpolate at all points in H* . On the other hand, X*f can be used 
to derive an truly interpolation function for points on the fundamental simplex, which will be 
developed in the next section. Recall 5k defined in (3.151. 

Theorem 3.20. For n > and f e C(U H ), define 



Kf(t) := £ /((^Kj,„(t), 



whe 



£j, n (t) = K(t and K(t) = jj-^- d £ 

1 ; keH; 



<£°<Mt). 



Then T*f € T n and it satisfies 



(3-23) Kf(xd^) = 



f( (d+l)n)' J ^ 



E%M' J £ H* \ H°. 
Furthermore, $* (t) is a rea/ function and it satisfies 

J- 1 *^ /CN d+i "•' 



w/tereN^ := N d+1 \ {j}. 
Proof. By definition, 

^™((d+l)n) = $ n((j+i)„) = ^ + X C l ^((d+ijn)' 



Since tiles R 1 ^ 1 , there exist m, 1 e Z^ 1 such that ( d ™)„ e and k-j = m+(d + l)nl. 
Thus, by Theorem [3TT7] 

= (<t>m, (f>o)^ = (<t>m,(f>a) = S mfi . 

Equivalently we can write the above equation as 

(3-25) = J 1 ' k =j + ( d+1 H leZ H +1 ' 

v ; J ' n M*H)»' wkivj/n otherwise; 



which proves (3.23). 
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To derive the compact formula for $* (t) we start with the following symmetry argument, 

i^E E Mt°) = 4 E E ^-(t) 4EE <Mt) 

= w\ E E *k(t) 4E E ^(*) + E **(t). 

Using the fact that Uo-eS ^o' J(T = an< ^ ( |3. 14| > , we deduce that 



E E Mt) - E E E «t) = E E i ^^( t ) 

<reSke(H„\H°)<r 0< i j<i+i<d+l<reekeHj l ' J cr 0<i,i«+J<«J+l kdHfr* 

= (d+i)i X! E 4"Vk(t) = (rf+i)! 53 4 n Vk(t), 

0<i,j<i+j"<d+l kgjj^i keH*\H° 

where the factor ^/v^'P 1 in the third summand comes from the fact that, by Lemma 



(i+j) 



J2aeg SkeH^<r summarizes over + x (.^^i^+i-i-j)! terms whereas S k eH*,' J summarizes 



3.8 



(<*+!)! 



terms. As a result, we obtain that 



i4 E E Mf) = E c k l Vk(t) + 53 k (t) = 53 4™Vk(t). 

' ' crGakGH„ keH*\H° keH° keH* 

In other words, we have shown that 

= w\ E with ^ = E «Mt<7). 

' ' tree keH„ 



We now evaluate the partial sum D n (t). Using Lemma 3.10 and the fact that t\+., -+td+i = 
0, we obtain 

D n (t) = J2 E M t )=E E e^S«i<i( fc *-**)*'+E,<,<-H- 1 (*«-*i)* I ) 
J =1 keHi j} i^iken^ 1 

d+lj-l n d+1 n-1 . ci+1 

_ TT gZirt^ TT ^ Znivtt = \ ^ TT sm7rnt ' c 7r»(n+l)tj TT ^^^ ^.(..-Dt, 
Z^ll^ 11 Z^ll sm7 rf 11 sinTrt; 

3=1 i=l l/=l i=j+l ^=0 3 = 1 i=l (=3 + 1 

■L-^- sin7rf; ^-^ sin7rnt. 



i=l 3=1 



''3 



Using symmetry and t\ + . . . + i^+i = again, we further derive that 

., 1 d+l . d+1 . 

^E 5 «(^-^n s ^E in ^r e " 7rmtj E ie/l><lWi e ^ +2E<elt<) - 

Since H* is symmetric under the mapping n : t i— > {t\, . . . , ij-i, — U, tj+i, . . . , td+i) and the 
definition of c k ™^ implies that c^. = , we see that the function X^keH* c k™Vk(t) is invariant 
under the action of Tj. As a consequence, we this function must be real valued. Hence, taking 
the real part in the above expression, we conclude that 

d+l 



E - E ( n 1^) ™^ E » s '< ( ' + 2 E'0. 

keH* 3=i v i<i<<i+i 1 / , rsl fi) v >' iei 

which completes the proof. □ 
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We note that the condition t\ + . . . + <d+i = implies various relations between trigonometric 
functions sin7rfi and cos7rfi. For example, in the case of d — 2, we have that 

sin 2t\ + sin 2t 2 + sin 2t% = —4 sin t\ sin ti sin t% 
cos 2ti + cos 2t 2 + cos 2t 3 — 4 cos ti cos i 2 cos t 3 — 1. 



Using these relations, it is possible to rewrite the compact formula in (3.241 in other forms. For 
instance, in the case of d= 2, the compact formula ( |3.24 ) becomes 

1 sin nnti sin -nnti 



371^ sin7r<i sin7rt2 
1 sin 7mi 2 sin imt^ 
3n d sin 7rf 2 sin irt^ 
1 sin 7m<3 sin nnti 



i nnts (| cos nt 3 + | cos n(ti — t 2 )) 
i nnti (| cos irti + | cos n(t 2 — £3)) 



+ CQS7rnt2 (2 COS7T f 2 _|_ l COS 7r(i 3 - ti)) , 
3n d sm7rf 3 sui7rti VJ J ' 

whereas a different formula is derived in |18j . given in terms of 0„(t). 

Using the compact formula of <&* , we can estimate the operator norm of I*. 

Theorem 3.21. Let ||/*||oo denote the operator norm of I* : C(Qh) >— > C(£Ih)- Then there is 
a constant c, independent of n, such that 

\\i:\\oo<c(logn) d . 

Proof. A standard procedure shows that 



^lloo = max ^ |*n(* ~ XdTV^) 
ten- 1 



Using the compact formula of $* in Theorem [3^20 it is easy to see that it suffices to prove that 



li '■= -7-, ; — j max 

3 (d+l)n d te u H 



e n 



k ... 1 :i< L d+i sin?r(ti (d+i) r 



< c(logn) d , 



1 < i < d + l, n > 0. 



prf+i 

if ' 



let := (ti, . . . , . . . ,*d+i) G R d . It is easy to see that {t^ : t G 



For te 

fi#} C [-1, Hence, enlarging the domain M* to {k G : -(d+l)n < fc; < (d+l)n, ki 

(mod d + 1), 1 < i < d}, we see that 



1 



n 



/, < 7- : — 7 max FT > 

J ~ (d + l)n d te[-i,i]* H , 

£—1 fc; — — n 



siii7m(t ; - 



< 



1 te[— 1,11 \ ^ 



E 

fc— — n 



sinn(ti 
sinnn(t — -) 



shi7r(t - -) 



< c(logn) d , 



where the last step follows from the usual estimate of one variable (cf. [34l Vol. II, p. 19]). □ 

We expect that the estimate is sharp, that is, ||/*|| > c(logn) d , but do not have a proof at 
this point. 



4. Discrete Fourier analysis on the simplex 
The fundamental domain of the lattice Ad is the union of the images of the fundamental 



simplex Ah under the group Q, as shown in (3.3 1. Hence, if we consider invariant functions 



under the group Q, then the discrete Fourier analysis on the fundamental domain in the previous 
section can be carried over to the analysis on the simplex Ah, which is developed below. 
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4.1. Generalized sine and cosine functions. In the case of one- variable, the invariant and 
anti-invariant sums of the exponential functions are cosine and sine functions. We now consider 
their analogous in our setting. It turns out that these functions have already appeared in 
the literature, as mentioned in the introduction. It should be pointed out, however, that our 
study is on the discrete Fourier analysis, which has little overlap with the previous study in 
the literature. Most of the overlap will appear in the next section, when these trigonometric 
functions are transformed to Chebyshev polynomials. 

Recall that the reflection group Ad is the permutation group Q. Denote the identity element 
in Q by 1. It is easy to see that 

2 



(4.1) 



o-ijCTjkVij = Oik, hj,k G N d+1 . 



For a G G, let \a\ denote the number of inversions in a . The group Q is naturally divided 
into two parts, Q + := {a G Q : |<r| = (mod 2)} of elements with even inversions, and Q~ := 
{(J E G ■ |cr| = 1 (mod 2)} of elements with odd inversions. The action of a G Q on the function 
/ : R^ 1 i — > R is defined by a fit) :— f(ta). A function / in homogeneous coordinates is called 
invariant under Q if af — f for all a G G, and it is called anti- invariant under Q if af = p(a)f 
with p(a) — 1 if a G G + and p(a) = — 1 if a G G~ . The following proposition follows immediately 
from the definition. 

Proposition 4.1. Define two operator V + and V~ acting on f(t) by 



(4.2) 



V ± f(t) := 



1 

(d+i)i 



]T f(ta) ± £ f(ta) 

aeg+ <?eg- 



Then the operators V + and V are projections from the class of H -periodic functions onto the 
class of invariant, and respectively anti-invariant functions. 



Recall that </>k(t) = e d + 1 Applying the operators V^ 1 to these exponential functions gives 
the basic invariant and anti-invariant functions. 



Definition 4.2. For k G H define 
TC k (t) := 7>+0k(t) 

TS k (t) := V'M*) = 



1 



(d+1)! 
1 



ik(t£T)+ ^ k ( tCT ) 

aeg- 

) k (ttr) - ^ </>k(tcr) 



aeg- 



and call them generalized cosine and generalized sine, respectively. 

By definition, TCk is invariant and TSk is anti-invariant. Because of the symmetry, we only 
need to consider them on the fundamental simplex A# defined in (3.2 1 or any other simplex 
Ajjct, a G G, that makes up . We shall work with A# below and recall that 



A 



H 



{te 



4 +1 



< ti - tj < 1, 1 < i < j < d + 1} 



In the case of d = 2 and <i = 3, it is an equilateral triangle and a regular tetrahedron, respec- 
tively; these regions are depicted in Figure 4.1, in which the corners are given in homogeneous 
coordinates. 

By definition, 0k and 0ko- with a G G lead to the same TCk- In fact, we have 
(4.3) TC kff (t) = TCk(ter) = TC k (t) for t G A H and a G Q. 

Thus, when working with TCk, w e can restrict k to the index set 

A := {k G H : fci > fc 2 > ■ ■ • > k d+1 } . 
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Figure 4.1. Reference triangle and tetrahedron. 



As for TSk, it is easy to see that TSk CT (t) = TSk(ter) = TSk(t) for a G G + and TSk CT (t) = 
TSk(tcr) = — TSk(t) for a G Q~ . In particular, TSk(t) — whenever two or more components 
of k are equal. Thus, when working with TSk, we only need to consider k G A°, where 

A° := {k e H : fci > fc 2 > • • • > k d+1 } , 

which is the set of the interior points of A. To describe the points on the boundary of A, we need 
to consider the compositions of the integer d + 1. A composition p of d + 1 is a decomposition 
of d + 1 , such that 

p= (pi,...,p e ) Pi>0, l<i<£<d+l, \p\:=p 1 + ...+p e = d+l, 

where I = £{p) is the length of the composition. Notice that the order of pi matters, different 
orderings are deemed to be different compositions, which is the difference between a composition 
and a partition. We denote the collection of compositions of d+ 1 by Cd+i] that is, 

C d+ i := G Z t{p) : 1 < Pi < \p\ = d + 1 for 1 < i < £{p) < d + l} . 

For p e Cd+i we further define 

A" := {k = ({fcif 1 , {fc 2 f 2 , • • • , {k e } pe ) G M : ki > kj ,1 < i < j < £} . 

Then evidently A° = A« d+1 and A = [j peCd+1 A p - 

Recall that kCJ = {kcr : a e Q} is the orbit of k under Q. The definition of A implies that, 
for k, j G A, kg n \Q = whenever k ^ j. It follows that 



p / PiW----Pe l - 
Since is isomorphic to kQ, where is the stabilizer of k, we also have 

( 4 - 4 ) TCxW = 17* E E ^(t) = |4 E ^(*)- 

1 'jekefreSj 1 y| jekG 

The length f = £(p) of p G C<j+i determines how many indices in k G A p are repeated, 
which determines the dimension of the boundary elements of A. For instance, if d = 2 then 
g 3 = {(1,1,1), (1,2), (2,1), (3)} and 



|k<?| 




A 1 ' 1 - 1 , 
A 1,2 U A 2,1 , 
A 3 , 
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where A e and A" consist of points in A that are on the edges and vertices of A, respectively. If 
d = 3 then Q = {(1, 1, 1, 1), (1, 1, 2), (1, 2, 1), (2, 1, 1), (1, 3), (3, 1), (2, 2), (4)} and 



w = < 



24, 

12, 

6, 

4, 

1. 



ke A° := A 1 ' 1 ' 1 ' 1 , 

k e tJ := A 1 ' 1 ' 2 U A 1 ' 2 * 1 U A 2 * 1 ' 1 , 

k e A 6 - 1 := A 2 ' 2 , 

ke A 6 - 2 := A 1 ' 3 U A 3 ' 1 , 

k = e A v := A 4 , 



where A-^, A e and A v consist of points in A that are on the faces, edges and vertices of A, 
respectively. 

We define an inner product on A# by 

(d+1)! 



</,<?> 



A H ■= 



1 

|Ah| 



f(t)g(t)dt = 



/(t)fl(t)tft. 



A, 



If /<? is invariant under G, then it follows immediately that (/, g) — {f,g)A H , where (•, •) is the 
inner product defined in (3.6 1 over . The generalized cosine and sine functions are orthogonal 
with respect to this inner product. 

Theorem 4.3. For k, j G A. 

(4.5) (TC k ,TCj) Af 
and for k, j e A° . 

(4.6) (TSk, TSj) a 



Ok.l 

W\ 



■pv- 



(d + 1)! 



'k,i, 



ke A p , 



(d+1)! 



Proof. Both of these relations follow from the identity (/,<?) = (f,g)A H f° r invariant functions. 
For (4.5), the invariance is evident and we only have to use the orthogonality of 0k in (3.51 and 
(4.4). For (4.6), we use the fact that TSk(t)TSj(t) is invariant under Q and the orthogonality 
of 0k on fi^f. □ 

The definition and orthogonality of these trigonometric functions have appeared in the liter- 
ature, we refer to |3J and its extensive references. However, the study in the literature is more 
on the side of algebraic polynomials, as will be discussed in Section 5 below. 



4.2. Discrete inner product on the simplex. By Theorem 3.17 and ( |3.5[ ), {0k : k S H ra } 
is an orthonormal set with respect to the symmetric inner product (•, •)* . Using the symmetry 
and the invariance of TCk and TSk under G, we can deduce a discrete orthogonality for the 
generalized cosine and sine functions. We define 



(4.7) 



A„ := H* n A = {k e M : k d+1 < k d < . . . < k x < k d+1 + {d+ l)n} . 



For n = 4, the point sets for d — 2 and d = 3 are depicted in Figure 4.1. 

By the definition of H*, the set { ^+i)n : ^ e ^} con tains points inside Ah- We will 
also need to understand points on the boundary of Ah, which can come from two types of 
points in H* . One part of the boundary points in Ah are also the boundary points of H* , 
whereas another part of the boundary points in Ah are points inside H° but on the faces of 
Ah- Accordingly, for p S C d +i, we define 



k e + (d + l)n > ki > k 2 > ■ ■ ■ > h} , 
k e + (d + l)n — h > k 2 > ■ ■ ■ > k t } . 
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Figure 4.2. A„ for d = 2 and d = 3 



Evidently, A° = A°' {1} and 

A n = |J (A°*UA»*). 

peCd+i 

If P 7^ then A° ,p is a [l — l)-face of A„ which, however, is a subset of H°, whereas A^ ,p 

is a (£ — 2)-face of A„ which is also a face of H* . More precisely, considering the orbits of the 
points in A„, we see that 

(4.8) H° = |J A° n vg and H* \ H° .= |J A^£. 

peCd+i pGC<i + i 

For d = 2, the sets of interior points, edge points and vertices of Ah are given, in this order, 
explicitly by 

A"' 1,1 ' 1 ={kei: fc 3 + 3n > k x > k 2 > k 3 } , 

A°-^ 2 U A"' 2 ' 1 U A^' 1 ' 1 ' 1 = {(2fc, -k, -k), (k, k, -2k), (3n - k,2k - 3n, -fc) : < fc < n} , 
A°< 3 U A^< 2 U A^ 2 - 1 = {(0, 0, 0), (2n, -n, -n), (n, n, -2n)} . 

In this case, the geometry in Figure 4.3 shows clearly that two of the three edges, A"' 1 ' 2 U A°' 2,1 , 
of the triangle Ah are not edges of the hexagon, and the vertex A°' 3 = {(0,0,0)} is not an 
vertex of the hexagon. 




Figure 4.3. Reference triangle Ah inside the hexagon 
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For d — 3, the boundary sets are given in detail in [TO] and they are denoted by A°, A£, 
A^' 1 , A^' 2 and A£, which correspond to points of interior, faces, two type of edges, and vertices. 
In terms of the sets A°' p and A^' p , they are expressed by 

Ki = A"' 1 ' 1 ' 1,1 , 

Af n = A"' 1 ' 1 ' 2 U A" 4,2,1 U A"' 2 ' 1 ' 1 U A^' 1 ' 1 ' 1 ' 1 , 
A^^A^U^UA^UAW, 

a; = a°' 4 ua^ 3 ua^ 2 ua^ 3 ' 1 . 

We now define a discrete inner product (•, -)A,n by 

1 



(/)5>A,n = 



(d + l)n d .^V 



^ A j ( " ) /( (d+i)„)s , (( £ ;+i)„)- 



where, with cj n ' defined in Definition (3.161, 



(d + 1)! 



(4.9) 



(n) 



» 



d+1 



(d + 1)! 



j G A»*. 



(Pf +Pi)!p 2 ! 

Let us verify the second equal sign in ( |4.9| . The case j G A° ,p is easy, since then j G H° so that 
Cj = 1 and Aj = ( d+1 ) ■ In the case j G A^ ,p , we have ji—ji = (d+ l)n which implies ji —ji = 
(d + l)n for i G I Pl and / G J P€ , where J Pl = {1, . . . ,pi} and J P(! = {d + 1 — pe + 1, . . . , d + 1}. 
Consequently, , , , J n G by the definition of , which implies that j G BP 1 *' so that 



» 



l/( Pl + w ) and jOl. 



In the case of d = 2 and d = 3, the values of Aj are given by 



in) 




if d 



and 



(") 



if d 



24, jGA°, 
12, j G A£, 
6, j G A*' 1 , 
4, j G A e „< 2 , 
l, jeA", 

We denote by TC n and TS'n the spaces of the trigonometric polynomials 

TC n := span {TC k : k G A„} , TS n := span {TS k : k G A° } , 
respectively. Since A„ contains integer points in a regular simplex, it is easy to see that 



(4.10) dimTC„=|A n | = 

Theorem 4.4. For fg G TC 2n -i 



n + d 
d 



and dimTS*„ = |A°| 



71 — 1 

d 



(4.11) (f,g)A H = (f,g)A, n . 

Moreover, the following cubature formula is exact for all f G TC^n-Xt 



(4.12) 
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In particular, 
(4.13) 



<TC k ,TC 



A,. 



, k,j e A r , 



Proof. Let / be a function invariant under Q. According to (|4.8[), 



d+l 



E 4 f((d+l)n) - E 

j£H° p£C d+1 
E C j ^((d+ljn) = X! 

Adding these two expressions together, we conclude that 



E 4 ^( (d+l)n)' 



d+l 



9,p 



(4-i4) E c i n) /(^)= E 

j6H; pec d+ i 



= E ^ 

jeA„ 

Replacing / by /g, we have proved that {f,g)* n = (/,.<?) a. n whenever fg is invariant. Hence, 
(4.111 follows from Theorem 13.17 Furthermore, since p^-j- /(t)dt = rgj^ J A ^ f(t)dt for 



d + 1 

P 

d+l 
P 



E C j f((d+l)n) + E C j 



jeA^ 



e,j> 



/((d+l)n) — E "\i -^((d+l)n)- 
JGA„ 



all invariant /, (4.121 follows from Theorem 3.18 

Furthermore, replacing / by TCkTCj in (4.141, we derive by (4.3 1 that 

<TC k ,T Cj ) A ,„ = (ZT 1 1) ^ E 4 B) TC k (^)TCj(^) 
1 1 



(d+ l)n d 


\Q\ 


1 


i 


(d+ l)n d 


\Q\ 


1 


\Q\ 


(d+ l)n d 


\Q\ 



E c i (n) E^((3T^) TC j((3TIF) 



leH* 



E c r } E^i^ra^) 



leH* 



tree 



IWlil 51 c i ( "Vk( T 2TiF)TC j ( 7 ^ IF ) = (0 k ,TCj);. 



leH* 



Using (3.25) and abbreviating k = j mod (d + 1)Z^ +1 as k = j, we further deduce that 



(TC k ,TCj) A ,ri = T^T E^ k '^ CT 



1 



(d + l)! 



(d+l)! 
{cr e g : k(7 = k} 



{a EG :j<J = k} 



A 



(n)> 



where the last equality follows readily from the definition of A„ and (4.9 1. 



□ 



Since the above proof applies to invariant /, it also applies to /, g G TS n since fg is invariant 
if both / and g are anti-invariant. Recall that TSk(^pj^) = when j S A n \ A° , we have also 
proven the following result. 

Theorem 4.5. Let the discrete inner product (-,-)a o ,?i be defined by 

d! 



(f,g)A",n - E /((d+l)n)f((d+l)n)- 
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Then 

(/,S>A> = (/,s}a h , f,geTS n . 

4.3. Interpolation on the simplex. Using invariance and the fact that the fundamental 
simplex Ah is the building block of the fundamental domain, we can also deduce results on 
interpolation on the fundamental simplex. For d — 3 the results in this subsection have appeared 
in [19j and the proof there can be followed verbatim when 3 is replaced by d. Thus, we shall 
omit the proof. Recall that the operator V ± is defined in (4.2 1. 

Theorem 4.6. For n > 0, and f £ C(A H ), define 

£ n /(t) := £ /( I3 J IF )^ n (t), £? n (t) := ^±11! £ TS k (t)TS k (^). 
jeA° keA° 

T/ien £„/ is the unique function in TS n that satisfies 

£n/((d+l)n) = /((d+l)n)' J G ^n' 

Furthermore, the fundamental interpolation function £° }n is real and satisfies 

S,*(t) = ^n" [e„(t - (3^) - e n _r(t - 

where means that the operator V~ is acting on the variable t and 9„ is defined in (3.181. 

The function C n f interpolates at the interior points of A„. We can also consider interpolation 
on A„ by working with the operator X*f in Theorem 3.20 which interpolates / on H° but 
interpolates a sum over congruent points on H* C H° . Using invariance, however, I*f leads to 
a genuine interpolation operator on A„ . 



Theorem 4.7. For n > and f G C(A#) define 



A< n > 



C/(t) := E /(ufc)^n(t), £„(t) := E AL n) TC k (t)TC k (^) 

T/ien £* / is i/ie unique function in TC n that satisfies 

£n/((d+l)n) = /((d+l)n)' J G 

Furthermore, the fundamental interpolation function i, n is given by 

# B (t)=Aj (B) ^j.n(t), 



where £^ n is defined in Theorem 3.20 and has a compact formula. 

Let ||£ n || and ||££|| denote the operator norms of £„ and £* , respectively, both as operators 
from C(Ah) >— > C(Ajj). From Theorems 4.6 and 4.7 an immediate application of Theorem 
|3.21| yields the following theorem. 

Theorem 4.8. There is a constant c independent of n, such that 

||£„|| < c(logn) d and ||£;|| < c(logn) d . 

It should be pointed out that the interpolation functions defined in these theorems are 
analogous of trigonometric polynomial interpolation on equally spaced points |34l Chapt. X]. 
These are trigonometric interpolation on equal spaced points in the simplex Ah, which can 
be easily transformed to interpolation on regular simplex, say {y : < yd < • ■ ■ < yi < 1} in 
R d . These interpolation functions are real, easily computable from their compact formulas, and 
have small Lebesgue constants. They should be the ideal tool for interpolation on the simplex 
in R d . 
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5. Generalized Chebyshev polynomials and their zeros 

The generalized sine and cosine functions can be used to define analogues of Chebyshev 
polynomials of the first and the second kind, respectively, which are algebraic orthogonal poly- 
nomials of d- variables, just as in the classical case of one variable. These polynomials have been 
defined in the literature, as noted in the Introduction. In the first subsection we define these 
polynomials and present their basic properties. Most of the results in this subsection are not 
new, however, we shall present a coherent and independent treatment, and some of the results 
on recurrence relations appear to be new. In the second subsection, we study the common zeros 
of these polynomials and use them to establish a family of Gaussian cubature formulas, which 
exist rarely. 

5.1. Generalized Chebyshev polynomials. The generalized trigonometric functions can be 
transformed into polynomials under a change of variables z : Wj^ 1 — > C d , defined by 

(5.1) «j=2i(t):=T(v(t) = n^^e^ J,t ) fc = 1,2, . 

j£v k Q 

where v fc = ({d+ 1 — fc} fe , {— k} d+1 ~ k ), and (k = 1, 2, . . . , d) are vertices of the fundamental 
triangle A# defined in Section 2. It is easy to see that Zk = Zd+i-k- The homogeneity of t 
shows that v fc • t = (d + l)(ti + . . . + tk) and, consequently, 

^ * (T)Jc 

\j\=k 

which shows that z%, . . . , Zd are the first d elementary symmetric polynomials of e 2mtl , . . . , e 27rltd + 1 
The same change of variables are used in [3] . 

Since TCk(t), k € A, is evidently a symmetric polynomial in e 21 ™* 1 , . . . , e 27 ™ td+1 , it is a 
polynomial in zi,...,Zd- Some of its properties can be derived from the recursive relations 
given below. 

Lemma 5.1. The generalized sine and cosine functions satisfy the recurrence relations, 
(5.3) TC J (t)TC k (t) = ^ T ^TC k+jCT (t) ) j,keA, 



(5-4) TCj(t)TS k (t) = i^TS k+j(T (t), j,keA, 

(5.5) TSj(t)TS k (t)=^^(-l)HTC k+jff (t), j,keA. 

Proof. From the definition of TC k , we obtain that 

Tc k (t)Tc k (t) = ±- J2 <Mt) x l- ^a(t) = tLJ2J2 ^r+j*(t) 

' y ' aeg ' y ' o-ee t eg a eg 

= I^i2 Yl 0(k+j<Ti-i)r(t) = T^ja ^(k+ja)r(t) 

11 reg^eg 11 t eg a eg 

= |CT 51 TC k+j CT (t), 

11 aeg 



which proves (5.3 1. The other two relations, (5.4) and (5.5 1, can be established similarly. □ 



The polynomials defined by TC k (t) under (5.1l are analogue of Chebyshev polynomials of 
the first kind, to be defined formerly below. We will also define Chebyshev polynomial of the 
second kind, for which we need the following lemma. 



;so 
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Lemma 5.2. Let v° := ^~ 2 f+ 1 ) , . .. , r^±% ^ en 



(5.6) TS v0 (t) = -^- n (e 2 ^-e— ) = II -^-^ 
Furthermore, 

kHd + 1 — kV 

(5.7) z fc TS v o(t)= TS v o +v ,(t), \<k<d. 
Proof. Let /3 := (d, d - 1, . . . , 1, 0) € Nf| +1 . Then 

e 3¥T v °- t = e -7rid(ii + ...+t<i + i) e 27ri/3.t _ e 27ri/3.t_ 

Using the well-known Vandermond determinant (see, for example, [211 p. 40]), 
P^)^^ 13 ) = det(x l d+1_J )i<. iJ < rf+ i = J| (xi-Xj), 



oes d+1 i<i<j<d+i 

3 



and setting Xj — e 2 ^ ltj , we obtain 

Tvw^E^^ n 

Furthermore, since t\ + . . . + i^+i = 0, the above equation immediately gives 
TS v o(t) = jj-^Yy II (e ftrt -"-e* rt -M) > 



Next we prove (5.7 1. Using (5.4), we have 

z fc TS v o(t) = TC v fc(t)TS v °(t) = 7^7j TS v o +vfc(T (t). 



from which the second equal sign of (5.6) follows 

Assume 1 := v fe a ^ v fc for some cr £ Q. By the definition of v fe , there exists an integer 
1 < i < d such that Zj = — fc and Zj + i = d + 1 — k. Consequently, v° + lj = ( rf + 2 ~ 2 J')( d + 1 ) —k = 
(d-2j)(d+i) + _ fc) = «° +x + Z i+1) which implies that v°+v k a E dA and TS v o +vfc(T (t) = 0. 
Thus, by the definition of the stabilizer Q v k, 



tc i ^ V* tc k\{d + l — k)\ c 

z fc TS v o(t)= p-pyyy 2 , T5 v o +v fc g (t) = — ^Tpryj — TS v o +vft (t), 

as the stabilizer of v fc has cardinality kl(d + 1 — k)l. □ 
The same argument that proves ( |5.6[ ) also shows that, for each k S A, 
TS k+v o(t) = dct ( x 4 Aj+/3 ) , A := (fci - k d+1 ,k 2 - k d+1 ,.. ■ ,k d - k d+1 ,0), 

where Xi — e 27 ™ 4i . The definition of A shows that A is a partition. Hence, according to [2TJ p. 
40], TSk+ v °(t) is divisible by TS V ° in the ring Z[x±, . . . ,x d ], and the quotient 

s\(xi, ...,x d ) = TS k +v°(t)/TS v o(t) 

is a symmetric polynomial in x\, . . . , x dl which is the Schur function in the variables x\, . . . ,x d 
corresponding to the partition A. Since this ratio is a symmetric polynomial, it is then a poly- 
nomial in the elementary symmetric polynomials z\, . . . , z d ; more precisely, it is a polynomial 
in z of degree (ki — k d+ i) / (d + 1) as shown by the formula [211 (3-5)]. This is our analogue of 
Chebyshev polynomials of the second kind. 

To simplify the notation, we find it convenient to change the index and define the Chebyshev 
polynomials of the first and the second kind formally as follows: 
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Definition 5.3. Define the index mapping a : A 



(5.8) 



a,(k) 



A', 



H+l 



1 



1 < z < d. 



Under the change of variables (5.1l and (5.8 1, define 



T a {z) := TCk(t) and U a {z) 



TSk+ v ° (t) 



a £ Ng. 



TS v o(t) 

FFe ca// T a (z) and U a (z) Chebyshev polynomials of the first and the second kind, respectively. 

It is easy to see that the mapping (5.8 1 is an isomorphism. Indeed, since ki + . . .+ fc^+i = 
for k G A, it is easy to see that the inverse of a is given by 

dp i— 1 

ki = a v - (d + 1) a v , l<i<d + l. 

We also note that a(v k ) — ■= ({0} x , 1, {0} d ~ k ), the fc-th element of the standard basis for 
the Euclidean space K d . 

In the recurrence relation that we shall state in a moment, we will need the definition of T a (z) 
in which a can have negative components. If a has a negative component, say on < 0, then 
ki < k i+ i, so that k in (5.8) does not belong to A. If k £ H, we define by k + the rearrangement 
of k such that k + £ A. By the definition of TCk, we have TCi(;+(t) = TCk(t) for all k £ H. 
Thus, if a has negative component, then we define 

T a (z) = T a +(z) where a + corresponds to k + by (5.8). 

Both T a (z) and U a (z) are polynomials of degree \a\ = a\ + . . . + ad in z. Moreover, both of 
them satisfy a simple recursive relation, which we summarize in the following theorem. 

Theorem 5.4. Let P a denote either T a or U a . Then 



(5.9) 



P a {z) = Pa 



and they satisfy the recursion relation 

'd+V 

i 



(5.10) 



ZiP a (z) 



X1 {z), a£Ni 



: P »+ a (j)(z), a£N d , 



in which the components o/a(j), j £ v l Q, have values in {—1,0, 1}, U a (z) = whenever a has 
a component o>i — —1, and 

1 < k < d 



T (z) = l, T £k (z) = z k , 



U (z) = l, U ek (z) 



d+1 

k 



1 < k < d. 



Proof. The relation (5.9 1 follows readily form the fact that — (ki — kj) = kj — ki and ( |5.8[ ). The 
relation (5.10) follows immediately from (5.3) and ( 5.4 1 . The values of Po and P tk follow from 
definitions and (5.7l. If a has a component a.i = —1 then, by (5.8l, ki(a) = k i+ i(a) — (d+ 1). 
A quick computation shows then ki(a) + v° — ki + i(a) + v° +1 , which implies that k + v° £ dA, 
so that TSk+ v °(t) = and U a (z) = 0. □ 

It is easy to see that |a(j)|, j £ v 4 , also takes value in { — 1, 0, 1}. As a result, in terms of the 
total degree |a| of P a , the right hand side of (5.10) contains only polynomials of degree n — l,n 
and n+ 1, so that it is a three-term relation in that sense. For d — 2, the relation (5.10) can 
be rewritten as a recursive relation, 

P a +et(z) = 3z 1 P a (z) - ^+(-1,1) (z) - P Q _ £2 (z), 

P a +e 2 (z) = 3z 2 P a (z) - P Q+(1 _i)(z) - Pc-aiz) 



:S2 
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which can then be used to generate a polynomial P a from lower degree polynomials recursively. 
The same, however, cannot be said for d > 3. For example, if d — 3 then (5.10 1 for k = 2 is 

6z 2 P a (z) = P a +e 2 (z) + P a Ml,-l,l)( Z )+ P a+(l,0-l)(z) + P a +(-lfi,l)(z)+P a -e 1 +P a +(-l,l-l)(z), 

which has two polynomials of |a| + 1 in the right hand side, P a + t2 (z) and P a +(i.-i,i)( z )- It 
is possible to combine the relations in (5.101 to write P a+Ei (z) = Ei=i (f^ 1 ) ZiP a (z) + Q(z), 
where Q(z) contains only linear combinations of {Pp} for — \a\ and |/3| = |a| — 1. In lower 
dimension, the exact forms of Q can be easily determined (for d = 3 see [29]), but the general 
formula for a generic d appears to be complicated and we shall nor pursuit it here. 

Next we show that T a and U a are orthogonal polynomials. The integral of the orthogonality 
is taken over the region that is the image of Ah with respect to a measure, or weight function, 
that comes from the Jacobian of the change of variables. Furthermore, since Zk i— > z<j+i-fc, we 
can consider real coordinates, denoted by x, 

_ z k + z d+l-k _ z k — z d+l-k 
x k — ~ j Xd+l-k — 7T. 

(5.H) 2 ! 21 

and x d+i = —j= z d+i , if d is odd. 
2 v 2 2 



for 1 < k < L|J , 



d+l 
H 



Combining (|5.1| and (|5.11|, our change variable becomes 

sin 2 7r(i M 



x. Let us define 



w(x) = w(x(t)) 



n 

l<H<v<d+l 



t v ). 



Lemma 5.5. 

(5.12) 



The Jacobian of the changing variable t t— > x is given by 

d 



det 



d(x!,x 2 , ...,x d ) 
d(t 1 ,t 2 ,. ■ ■ ,t d ) 



= 2 



n 



( d+i ) 

fe=l V k ) 



,(x)y- 



Proof. Let di denote the partial derivative with respect to t t for 1 < i < d. We first prove that 

d 



(5.13) 



det 



d{zi,z 2 , ■ ■■ ,z d ) 
d(t 1 ,t 2 , ■■■,t d ) 



2m 
fc=i V fc / 



n 



x ]J[ IsiiiTT^ - 1„)|. 
l</i<i/<d+l 



Using the Jacobian of the elementary symmetric polynomials with respect to its variables and 
t\ + . . . + td+i = 0, this can be shown as in (3j (5.9)]. We give an inductive proof below. 
Regarding t\, t 2 , . . . , t d+ i as independent variables, one sees that 

dzk 2ni 



E 



( d+1 ) 



For each fixed j, split I C N d +i as two parts, one contains {j, d+l} and one does not, so that 
after canceling the common factor, we obtain 



dzk 



dz k 
dt d+1 



Hence, setting fj k := 
we have shown that 



2ttj 

(T) 

2ni , 

£ 



E 

je/c N ^+ 1} ,|/|=fc 



E 

d+le/CN^ } ls |/|=fe 



e 27r *E„ e / ^ 



E 



|/|=fc-l 



7CN 



ay 



det 



0(zi,22, 



|7|=fe- 



1 e 2wl ^"e/ ** and defining the matrix F d := (ff k ) 1 



l,kJi<j,k<d' 



d(t 1 ,t 2 , ■■■,t d ) 



n ?Sy x n ( e2 "* j - e2 "* d+i ) x det (^)- 
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By definition, ff x = 1; splitting according to if the part contains d or not, it follows that 



f*k~flk= (e^'-e 2 **') X 



/CN^' d} ,|/| = fc-2 



which allows us to use induction and show that 

det(F rf ) = Yl (e 2 "*" ~ e 2 ™*") x dct ( F <i-i) = • • • = JJ (e 27 ™*" - e 27 "*" 

l<jtl<d— 1 l<[i<v<d 



This proves (5.13) upon using (5.6 K Now, a quick computation using (5.11) shows that 

d(xi,x 2 , ...,Xd) 



det 



d(x 1 ,x 2 , ...,x d ) 



d(ti,t 2 , ...,t d ) 



det 



d(z 1 ,z 2 , ...,Zd) 



det 



d(zi,z 2 , ...,z d ) 
d(x 1 ,x 2 , ■■■,x d ) 



fe=l V k ) l<n<v<d+l 



which is exactly (5.12) 



□ 



Under this change of variables, the domain A# is mapped to 

A* := I x = x(t) € R d : t e ]J[ suitt^ - i/) > 1, 

^ l<i<j"<d+l ' 

which is the image of A# under the mapping t i— > x. In the case of d = 2 and d — 3, A* is 
depicted in the Figure 5.1. 



M. 775.°) 





Figure 5.1. The region A* for d — 2 and d = 3 
We will need the cases of /i = — 1/2 and = 1/2 of the weighted inner product 

(f,9)w-=c a f{z)g{z)w a {z)dx. 



where c a is a normalization constant, c a := 1/ J. m w a (z)dx. The change of variables t 
shows immediately that 

r— r / /(t)cft = c_i / 

l A #l J A h 2 J A' 

where 

d fd+l\ 



1 



w 2 (z) 



A 



H 



det ( X 1 '^""'^ ) 

V d(t 1 ,t 2 ,...,t d ) I 



fc=i 
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Furthermore, by (5.6 1 and (15.5 



d{d+\) d(d+l) 



(*) = ("?) 2 [(d+l)!] 2 [TS v o(t)] 2 = (-|) ' (d+l)!X)(-l) H TC v . +v . <r (t) I 



which shows that w is a polynomial in z with a total degree of 2d, and furthermore, by (4.6 1 

2^2 d jd+i 



Ci 



1/ ^ ^ ^(t))^ = 1/ jT W (x(t))dt = y ^ n 



The orthogonality of T K and [/ K , respectively, then follows from the orthogonality of TCk and 
TSk with the change of variables. More precisely, Theorem |4.3| leads to the following theorem: 



Theorem 5.6. The polynomials T a and U a are orthogonal polynomials with respect to w 1 I 2 
and w 1 ' 2 , respectively, and 

(5.14) {T a ,Tp) w _i = d a 6 aif3 , d a = * if a € A p , 

v p J 

(5.15) (U a ,Ufi) wh =jh^r v «,^G< 



where A p is the image of A p under the index mapping (5.8), 

A p ={ae<:a = a(k), keA"}. 

The orthogonality of these polynomials has been established in the literature; we refer to j3j 
and the references therein. By (4.10), the set {T a : a £ Nq} and, respectively, {U a : a € Nq} 



forms a mutually orthogonal basis for the space H d of polynomials in o?-variables with respect 
to u> - 2 and, respectively id'. 

5.2. Zeros of Chebyshev polynomials of the second kind and Gaussian cubature. It 

is well known that zeros of orthogonal polynomials of one variable are nodes of the Gaussian 
quadrature. We consider its extension in several variables in this subsection. First we review 
the background in several variables. 

Let w be a nonnegative weight function defined on a compact set £1 in R d . A cubature 
formula of degree 2n — 1 for the integral with respect to to is a sum of point evaluations that 
satisfies 

N 

f(x)w(x)dx = Aj/(y ), Xj G R 
! j=i 

for every / G n d , n _ 1 . The points x J = {x\, x\, . . . , x J d ) are called nodes of and the numbers Xj 
are called weights of the cubature. It is well-known that a cubature formula of degree 2n — 1 
exists only if N > dimll^^. A cubature that attains such a lower bound is called a Gaussian 
cubature. 

In one variable, a Gaussian quadrature always exists and its nodes are zeros of orthogonal 
polynomials. For several variables, however, Gaussian cubatures exist only rarely. In fact, the 
first family of weight functions that admit Gaussian cubature were discovered only relatively 
recent in [5]. When they do exist, however, their nodes are common zeros of orthogonal poly- 
nomials. To be more precise, let P n := {P K : n € Nq, \k\ — n} denote a basis of orthonormal 
polynomials of degree n with respect to w(x)dx. Then, it is known that a Gaussian cuba- 
ture exists if and only if P„ (every element in P„) has dimll^^ real common zeros; in other 
words, the polynomial ideal generated by {P a : \a\ — n} has a zero dimensional variety of size 
dimll^^. See [TT], [23l [26] for these results and further discussions. 

Below we shall show that a Gaussian cubature exists for wi, which makes the second 
family of examples that admits Gaussian cubature. First we study the zeros of the Chebyshev 
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polynomials of the second kind. Because of Zk = Zd+i-k and (5.9 1, the real and complex parts 
of the polynomials in {U a (z) : \a\ — n} form a real basis for the space of orthogonal polynomials 
of degree n. Hence, we can work with the zeros of the complex polynomials U a (z). 

Let Y n and Y° be the image of j ^ d ^ n '■ j € A„ | and | ^ 1 ^ )n '■ j £ A° j under the mapping 
t i — ► x respectively, 

Yn^^ddiFn)-- jeA„} and Y° :— {z(^jy^) : j £ A°|. 

Theorem 5.7. The set Y°, d is the variety of the polynomial ideal (U a (x) : \a\ = n). 

n, then 

(d + l)n. Hence, in order to show that U a 



Proof. The definition of U a and (5.8 1 shows that if U a has degree n, that is, \a 
k e A° in TS v o + k/TS v o satisfies k\ — k d+ \ 



vanishes on Y° +d , it suffices to show, by the definition of v°, that 

TS k (pfi^) =0 for k £ A° and k x - k d+1 = (d+ l)n. 

Recall that a^j denotes the transposition that interchanges i and j. A simply computation 
shows that 



(d+l)n 



(d+l)n 



2i sin 



- jd+l)(h - fcd+l) 



e (d+iy 



-({ii+id+i){ki+k d+1 )+2Y.i = -2 3v k u) 



(d+l) 2 n 

which is whenever j £ A„ and k £ A„ satisfies k\ — k d+ i = (d + l)n. Hence, it follows that 



TSk ((d+x)n) - 2\g\ ^ £<T (( rf+Jl )™) ^ k ( J (rf+i|«) 



o. 



where e CT takes value of either —1 or 1. 

B y gg , = |A° +d | = 



□ 



5.7 



dimn^-p According to the general theory, Thco- 

implies that the Gaussian cubature exists for the weight function u>5 on A*. In fact, 
the precise form of this cubature is known. 

Theorem 5.8. For on A* , a Gaussian cubature formula exists, which is given explicitly by 
f(x)w 2 (x)dx = 

I A' 



(5.16) 



where the normalization constant ci = 



- rr 



Proof. The change of variables t i— > x in ( |5.l[ ) and (5.11 ) shows 

[(d+1)!] 2 



f(z)w 2 (z)dx = 



|Ah| 

dl(d + l)I 
(n + ef) d 



/(z(t))|TS v o(t)| z dt 

I . 2 

| TSv °((n+d)(d+l)) f( Z ( 



(n+d)(d+l) 



it)). 



where the last step follows from the fact that TS v °(t) vanishes on the boundary of A. This is 
exactly the cubature (4.121 applied to the function /(z(t)) |TS v °(t)| . □ 

The set of nodes of a Gaussian cubature is poised for polynomial interpolation, that is, there 
is a unique polynomial P in H^-i such that P(£) = /(£) for all £ £ Y° +d and a generic function 
/. In fact, let K n (-, •) be the reproducing kernel of in L 2 L , which can be written as 

W 2 

K n (x,y):= ^2{d+l)lU a (z)U a (w) 

\cx\—n 
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by (5.151, where x,y is related to z, w according to (5.111, respectively. Then it follows from 



the general theory (cf. [TT]) that the following result holds: 

Proposition 5.9. The unique interpolation polynomial of degree n on Y° +d in n„_j is given 
by 

(5.17) C n f(x(t))= f( X ( ( n +d)(d+l) )) W ( X ( (n+d)(d+l) )) K n(t, (n+d)(d+l) )» 



where K*(t,s) — K n (x,y) with x,y related to t,s by (5.1l and (5.11), respectively 



In fact, this interpolation operator is exactly the one in Theorem |4.6| under the change of 
variables t i— > x. In particular, we can derive a compact formula for K n (x,y) as indicated in 
that theorem. 



5.3. Cubature formula and Chebyshev polynomials of the first kind. For the Cheby- 
shev polynomials of the first kind, the polynomials in {T a : \a\ = n} do not have enough 
common zeros in general. In fact, as shown in [13], in the case of d — 2, the three orthogonal 
polynomials of degree 2 has no common zero at all. As a consequence, there is no Gaussian 
cubature for the weigh function w~ 2 on A* in general. 

However, the change of variables t i— * x shows that (4.12 1 leads to a cubature of degree 2n — 1 
with respect to ro~5 based on the nodes of Y n . 



Theorem 5.10. For the weight function w ? on A* the cubature formula 

' J2 X i n) f( z (^d+T)))-' v/ 6 n 2n _ l5 



(5.18) 



f(x)w 2 (x)dx = —r 

A* (d- 



Y)n° 



holds, where c_ 



Vrf+irfinLi and A j ,l) is 9 wen b y (El 1 



It follows from (4.10 1 that the \Y n \ = ( n ^" ) = dim 11^ and Y n includes points on the boundary 



of A*, hence, the cubature in (5.181 is an analog of the Gauss-Lobatto type cubature for w 2 
on u> - 5. The number of nodes of this cubature is more than the lower bound of dimll^^. 
Such a formula can be characterized by the polynomial ideal that has the set of the nodes as its 



variety. Indeed, according to a theorem in [33], the set of nodes Y n of the formula (5.181 must 
be the variety of a polynomial ideal generated by dimll^ +1 linearly independent polynomials of 
degree n+l, and these polynomials are necessarily quasi-orthogonal in the sense that they are 
orthogonal to all polynomials of degree n — 2. Giving the complicated relation between ideals 
and varieties, it is of interesting to identify the polynomial idea that generates the cubature 



(5.181. This is given in the following theorem. 



Theorem 5.11. Let v* = kT 1 ^*) = (d+ l,{0} d "\ -d- 1). For a = a(k) as in ([KsJ and let 

a* :— a ((k — v*) + ) . Then Y n is the variety of the polynomial ideal 



(5.19) 



(T a (x) -T( a ,)(x) : \a\ = n+l). 



Furthermore, the polynomial T a (x) — T/ a *\(x) are of degree n+l and orthogonal to all polyno- 
mials in II Tl _2 with respect to . 



Proof. From the definition of TCk, we have 

TC k (t)-TC k ^(t)^^E( 0k + t CTld+1 



9 keri -(i+ i+v* 



tree 



(ta). 
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A direct computation shows that 

</>k(t) + <Pka 1:d+1 (t) - 0k-v* (t) - 0far M+1 +v* (t) 

= e aTr(^^=2 2 ' £ - t -+( fe i+ fe <i+i)(*i+ t £ i+i)) ^ e T»(*i-*d+i) _ e T»(*d+i-*i)^ 
x ^ e ari( fc i _ ' c d+i _ ' i ~ 1 )(*i- t £i+i) _ e ari( fc i-' i -'£i+i-<i-i)(t < i+i-ti)~j 
= -4sin 7 r(i 1 - t d+1 ) shirr ( fcl ^ +1 - l) (i x - t d+1 ) e^&= ^^^wlfW^)) . 
Hence, for any j £ A n and k G EI with ki — fcd+i — (n + 1), 

(0k + 4>ka 1}d+1 - 0k-V - ^ko-i.d+l+v*) ( (d+l)„ 

which yields that 

(TC k - TC k _ v .) (p^) - for k G H with k x - k d+1 = (n+ l)(d+ 1). 

As in the proof of Theorem 5.7 by (5.8 1, this shows that T a — T a * vanishes on Y n . 

Moreover, suppose k £ A and set j := (k — v*) + ; since k% > ... > kd+i and fej = kj 
mod d + 1 for i > j, it follows that ji = fci — (d + 1) if fci > /c2 and Ji = fci if fci = 
fe 2 , and = fc d+ i + (d+ 1) if fc d+ i < fc rf and = k d+1 if = fc rf . Consequently, 
(ji —jd+i) — (ki — kd+i) £ {0, — d — 1, —2d — 2}, which shows that |a + | £ {\a\, \a\ — 1, |a| — 2}, 
so that T a * is a Chebyshev polynomial of degree at least n — 1 and T a — T Q * is orthogonal to 
all polynomials in n^_ 2 - □ 

We note that the general theory on cubature formulas in view of ideals and varieties also 
shows that there is a unique interpolation polynomial in 11^ based on the points in Y n . This 
interpolation polynomials, however, is exact the interpolation trigonometric polynomial in The- 
orem |4.7| We shall not stay the result formally. 



References 

[1] H. Bacry, Generalized Chebyshev polynomials and characters of GL(N, C) and 5L(N, C), Group Theoret- 
ical Methods in Physics, Lecture Notes in Phys., vol. 201, p. 483-485, Springer- Verlag, Berlin. 1984. 

[2] H. Bacry, Zeros of polynomials and generalized Chebyshev polynomials. Group theoretical methods in 
physics, Vol. I (Yurmala, 1985), 481-494, VNU Sci. Press, Utrecht, 1986. 

[3] R. J. Beerends, Chebyshev polynomials in several variables and the radial part of the Laplace-Beltrami 
operator, Trans. Amer. Math. Soc, 328 (1991), 779-814. 

[4] R. J. Beerends and E. M. Opdam, Certain hypergeometric series related to the root system BC, Trans. 
Amer. Math. Soc. 339, (1993), 581-609. 

[5] H. Berens, H. Schmid and Yuan Xu, Multivariate Gaussian cubature formula, Arch. Math. 64 (1995), 26-32. 

[6] J. H. Conway and N. J. A. Sloane, Sphere Packings, Lattices and Groups, 3rd ed. Springer, New York, 
1999. 

[7] A. Debiard, Polynomes de Tchebvchev et de Jacobi duns un espace euclidien de dimensron p, C. R. .Acad. 

Sci. Paris 296 (1983), 529-532. 
[8] K. B. Dunn and R. Lidl, Multi-dimensional generalizations of the Chebyshev polynomials, Proc. Japan 

Acad, 56 (1980), 154-165. 

[9] K. B. Dunn and R. Lidl, Generalizations of the classical Chebyshev polynomials to polynomials in two 
variables, Czechoslovak Math. J. 32 (1982), 516-528. 

[10] D. E. Dudgeon and R. M. Mersereau, Multidimensional Digital Signal Processing, Prentice-Hall Inc, En- 
glewood Cliffs, New Jersey, 1984. 

[11] C. F. Dunkl and Yuan Xu, Orthogonal polynomials of several variables, Encyclopedia of Mathematics and 
its Applications, vol. 81, Cambridge Univ. Press, 2001. 

[12] R. Eier and R. Lidl, A class of orthogonal polynomials in k variables, Math. Ann. 260 (1982), 93-99. 

[13] B. Fuglede, Commuting self-adjoint partial differential operators and a group theoretic problem, J. Func- 
tional Anal. 16 (1974), 101-121. 

[14] J. R. Higgins, Sampling theory in Fourier and Signal Analysis, Foundations, Oxford Science Publications, 
New York, 1996. 



38 



HUIYUAN LI AND YUAN XU 



[15] T. Koornwinder, Orthogonal polynomials in two variables which are eigcnfunctions of two algebraically 

independent partial differential operators, I - IV, Nederl. Acad. Wetensch. Proc. Ser. A77 = Indag. Math. 

36 (1974), 48-66 and 357-381. 
[16] A. U. Klimyk and J. Patera, (Anti)symmctric multivariate exponential functions and corresponding Fourier 

transforms. J. Phys. A 40 (2007), 10473-10489. 
[17] A. U. Klimyk and J. Patera, (Anti) symmetric multivariate trigonometric functions and corresponding 

Fourier transforms. J. Math. Phys. 48 (2007), no. 9, 093504, 24 pp. 
[18] H. Li, J. Sun and Yuan Xu, Discrete Fourier analysis, cubature and interpolation on a hexagon and a 

triangle, SIAM J. Numer. Anal, 46 (2008) 1653-1681. 
[19] H. Li and Yuan Xu, Discrete Fourier analysis on a dodecahedron and a tetrahedron, Mathematics of 

Computation, to appear (preprint, 2008). 
[20] R. Lidl, Tschcbyschcff polynomc in mehreren variabelen, J. Reine .Angew. Math. 273 (1975), 178-198. 
[21] I. G. Macdonald, Symmetric functions and Hall polynomials, 2nd ed, Oxford Univ. Press Inc., New York, 

1995. 

[22] R. J. Marks II, Introduction to Shannon Sampling and Interpolation Theory, Springer- Verlag, New York, 
1991. 

[23] I. P. Mysoviskikh, Interpolatory cubature formulas, Nauka, Moscow, 1981. 

[24] A. N. Podkorytov, Order of growth of Lcbcsguc constants of Fourier sums over polyhedra (Russian. English 
summary), Vestnik Leningrad. Univ. Mat. Mekh. Astronom, 1982, 110—111, 127. 

[25] P. E. Ricci, Chcbyshcv polynomials in several variables, (Italian), Rend. Mat. 11 (1978), 295-327. 

[26] A. Stroud, Approximate calculation of multiple integrals, Prentice-Hall, Englcwood Cliffs, NJ, 1971. 

[27] J. Sun, Multivariate Fourier series over a class of non tensor-product partition domains, J. Comput. Math. 
21 (2003), 53-62. 

[28] J. Sun, Multivariate Fourier transform methods over simplex and super-simplex domains, J. Comput. Math. 
24 (2006), 305-322. 

[29] J. Sun, A new class of three-variable orthogonal polynomials and their recurrences relations, Science in 

China, Series A, 51 (2008), 1071-1092. 
[30] X. Sun, Approximation on the Voronoi cells of the A^ lattice, manuscript, 2008. 

[31] V. A. Yudin, Lower bound of the Lebesgue constants. (Russian) Math. Notes 25 (1979), no. 1-2, 63 - 65. 
[32] Yuan Xu, On orthogonal polynomials in several variables, Special Functions, q-series and Related Topics, 

The Fields Institute for Research in Mathematical Sciences, Communications Series, Volume 14, 1997, p. 

247-270. 

[33] Yuan Xu, Polynomial interpolation in several variables, cubature formulae, and ideals, Advances in Comp. 

Math., 12 (2000), 363-376. 
[34] A. Zygmund, Trigonometric series, Cambridge Univ. Press, Cambridge, 1959. 

Institute of Software, Chinese Academy of Sciences, Beijing 100080, China 
E-mail address: hynliOmail.rdcps.ac.cn 

Department of Mathematics, University of Oregon, Eugene, Oregon 97403-1222. 
E-mail address: yuanSmath.uoregon.edu 



